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The First VLBI Detection of an Ultracool Dwarf; Implications for the Detectability of Sub-stellar 
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ABSTRACT 

We present milliarcsecond-resolution radio very long baseline interferometry (VLBI) observations of the 
ultracool dwarfs TVLM 5 1 3-46546 (M8.5) and 2MASS J00361617+1821 104 (L3.5) in an attempt to detect 
sub-stellar companions via direct imaging or reflex motion. Both objects are known radio emitters with strong 
evidence for periodic emission on timescales of about 2 and 3 hours, respectively. Using the inner seven 
VLBA antennas, we detect unresolved emission from TVLM 5 13^6546 on a scale of 2.5 mas 50 stellar 
radii), leading to a direct limit on the radio emission brightness temperature of Tg > 4 x 10^ K. However, 
with the higher spatial resolution afforded by the full VLBA we find that the source appears to be marginally 
and asymmetrically resolved at a low S/N ratio, possibly indicating that TVLM 5 1 3^6546 is a binary with 
a projected separation of 1 mas (~ 20 stellar radii). Using the 7-hr baseline of our observation we find 
no astrometric shift in the position of TVLM 5 13-46546, with a 3fT limit of about 0.6 mas. This is about 3 
times larger than expected for an equal mass companion with a few-hour orbital period. Future monitoring 
of its position on a range of timescales will provide the required astrometric sensitivity to detect a planetary 
companion with a mass of ^ 10 My in a > 15 d (> 0.06 AU) orbit, or with a mass of ^ 2 Mj in an orbit of 
> 0.5 yr (> 0.3 AU). 

Subject headings: stars: low-mass, brown dwarfs - radio continuum: stars 



1. INTRODUCTION 

In recent years unexpectedly strong radio emission has been 
detected from very low mass stars and brown dwarfs (here- 
after, ultracool dwarfs), revealing that these objects are capa- 
ble o f generating and dissipating kG-strength magnetic field s 
(e.g., BergeretiOJIZOOlt lBergei]l2006t iHallinan et al]l2008l) . 
The radio luminosity remains nearly uniform from early-M 
to mid-L dwarfs, even though other magnetic activity indica- 
tors (Ha and X-rays) decline by about two orders of magni- 
tude over th e same spectral type range Ce.g. JWest etaPlBol 
|Berge/2006). Thus, radio observations provide a particularly 
powerful probe of the magnetic field properties. 

The radio emission from several ultracool dwarfs also ex- 
hibits a periodicity that matches the stellar rotation veloc- 
ity (vsin/), thereby pointing to a large-scale, axisymmet- 
ric, and stable magnetic fi eld configuration on timescales of 
hours to years (iBereer et al. 2005; Hallinan et al. 2006, 2007; 
iBerger et al.ll2009l) These conclusions are also supported by 
observations of Zeeman broadening in FeH molecular lines 
(iReiners & Basri 2007) and t ime-resolved optic al spectropo- 
larimetry (" Donati et al.ii2008l:lMorin et al.ll2008l) . 

Since ultracool dwarfs are fully convective, the solar-type 
ail dynamo, which operates at the sheari ng interfac e be- 
tween the radiative and convective zones (P arke3fT955l) . can- 
not be responsible for generating and amplifying the inferred 
fields. However, recent numerical simulations suggest that 
large scale axisymmetric fields can still be generated in fully 
convective objects, at least for conditio ns that rou ghly corre- 
spond to mid-M dwarfs (M ~ 0.3 M0; ,Browning|2OO8) . Sim- 
ulations that correspond to the conditions in ultracool dwarfs 
are challenging and have not been investigated so far. As a re- 
sult, observational constraints on the scale, geometry, and ori- 
gin of the fields are essential. The spectroscopic Zeeman tech- 
niques are currently of limited utility for the dim and gener- 
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ally rapidly-rotating ultracool dwarfs, since they require very 
high signal-to-noise ratios and have reduced s ensitivity at high 
rotation velocities due to line broadening ( Reiners & Basril 
2007). Detailed radio observations are thus crucial. 

Here we present the first very long baseline interferometry 
(VLBI) observations of ultracool dwarfs, aimed at providing 
further constraints on their magnetic field properties. In par- 
ticular, these observations have sufficient angular resolution 
to detect astrometric shifts of a few tenths of a milliarcsecond 
that may be exerted by a sub-stellar companion with a period 
of a few hours. They also allow us to directly image a radio- 
emitting companion down to milliarcsecond scales. The pres- 
ence of such putative companions may enhance the magnetic 
activity through direct or tidal interaction. We detect the well- 
studied M8.5 dwarf TVLM513-46546 with the VLBA, the 
first such detection for any ultracool dwarf. The detectability 
of this object with the VLBA suggests that long-term astro- 
metric monitoring could reveal the presence of a planetary- 
mass companion on a > 15 d orbit. Such a detection would 
also usher in a new era of extrasolar planet studies through 
radio astrometry. 



2. TARGETS AND OBSERVATIONS 

The M8.5 dwarf TVLM 5 13^6546 (hereaf ter TV LM5 1 3- 
46) is located at a distance of 10.6 pc (Dahn et al. 2002n. 
It was first detected at radio wavelengths by .Bergen (i2002l) 
who found persistent emission and a circularly po larized flare 
lasting about 15 min at 8.5 GHz. Subsequently. lOsten et al.l 
(2006) detected quiescent radio emission at 4-.9 and 8.5 GHz, 
but did not detect any flares. iHallinan et al.l (12006!) found ro- 
tational modulation in the ra dio emission w i th a p eriodicity 
of ^ 2 hr, and subsequently, IHallinan et al.l (l2007l) detected 
periodic bursts of circularly polarized radio emission with a 
period of 1.96 hr. Most recently, Bergeretal. (2008) con- 
firmed the presence of polarized flares combined with quies- 
cent emission in simultaneous multi-wavelength observations 
(radio. X-ray, and Ha). These authors also found weak X-ray 
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emission and periodic Ha emission. Thus, the current picture 
of the radio emission from TVLM513^6 is one involving 
different activity states; quiescent radio emission appears to 
always be present (200-400 /iJy at 8.5 GHz), but variability 
and periodic bursts have not been detected in all observations. 

We also observed 2MASS J00361617+1821104 (hereafter, 
2M 0036+18), an L3 .5 dwarf located at a distance of 8.8 
pc (iDahn et al.ll2002i) . This obje ct has also b e en pre viously 
detected at radio wavelengths: iBerger et al.l (12005b found 
strongly variable and periodic radio emission with a period- 
icity of 3 hr, which was later confirmed by iHallinan et al.l 
([2008). The average 8.5 GHz flux density measured previ- 
ously was about 135 /iJy (iBerger et al.ll2005]) . 

2.1. Radio Observations 

TVLM5 13-46 and 2M0036+18 were observed simultane- 
ously with the Very Long Baseline Array (VLBA), the Very 
Large Array (VLA), and the Green Bank Telescope (GBT) on 
2008 March 30. Four hours of on-source time in the VLBI ob- 
servations were obtained at 05:45-12:32 UT for TVLM513- 
46 and 13:40-20:35 UT for 2M0036H-18; the VLA observa- 
tions covered the entire time range. All data were analyzed 
with the Astronomical Image Processing System (AIPS). 

2.1.1. VLA 

VLA observations were obtained at a frequency of 8.5 GHz 
with 26 antennas in the C configuration. We used the stan- 
dard continuum setup with 2 x 50 MHz bands and full po- 
larization. Flux calibration was obtained using J1331H-305, 
while J1455H-2131 (for TVLM513^6) and J0019H-2021 (for 
2M0036H-18) were used as complex gain calibrator sources. 
The bootstrapped flux densities for these two sources were 
0.081 Jy and 0.866 Jy, respectively. About once every 50 min 
additional scans of secondary cahbrators were obtained, re- 
sulting in corresponding gaps in the target source light curves. 

2.1.2. VLBA+VLA+GBT 

The VLBI observations included the 10-antenna VLBA, the 
GBT, and the VLA as a phased array. Phase-referenced ob- 
servations were carried out with a data rate of 256 Mbit s~' 
in dual polarization, using 2-bit sampling. Eight base-band 
channels of 8 MHz bandwidth each covered an aggregate 
bandwidth of 32 MHz. The correlator dump time was 2 s. 
The complex gain calibrators J1455H-2131 and J0019H-2021 
were located at distances of 1.8° and 4.4°, respectively. 

Unfortunately, the phase calibrators were not sufficiently 
bright to phase up the VLA between scans on the target 
sources, and the phased VLA data had to be removed from 
the VLBI experiment. Similarly, the GBT experienced bad 
weather, and the telescope had to be stowed for nearly 2.5 hr 
due to freezing rain during the second half of the experiment. 
As a result of the poor phases we also removed the GBT from 
the VLBI experiment. Finally, of the ten VLBA antennas, we 
removed Hancock (Hn), which experienced heavy snow fall 
and hence elevated system temperatures. We therefore sepa- 
rately analyze the VLA and VLBA data sets. 

3. RESULTS 

The VLA data allow us to determine the flux densities of 
our target sources, to construct light curves, and to constrain 
the search areas for counterparts on VLBI scales. 

3.1. VLA observations 



Both sources are detected in the VLA observations. 
TVLM513-46 has a flux density of 539 ± 19 ^Jy, and is 
located at a position of RA=15*'0 1^08. 162^ (±0.002 s), 
Dec=+22°50'01.673" (±0.035"). The synthesized beam size 
is 2.3" X 2.2". No circular polarization is detected in the in- 
tegrated data, with a 3(7 limit of < 9%. The light curve 
reveals that the emission is composed of a quiescent compo- 
nent and short-duration flares lasting for only a few minutes 
but with circular polarization approaching ~ 100% (Figure[T]). 
The most prominent flares are detected at 08:25 and 09:22 
UT. During the first two hours, when no flares are detected, 
the flux density is 445 ± 28 ^Jy. 

In light of previous results indicating that TVLM5 13-46 
exhibits periodic radio variability, we conduct ed a peri odicity 
search using the Lomb-Scargle periodogram (IScargldlT982n : 
Figure |2] We find a signal at a period of 58.8 min, exactly half 
of the previously reported period. Surprisingly, no obvious 
signal is detected at the harmonic of 1 17.6 min, even though 
three cycles at that periodicity fit into the observation. The 
58.8-min period is the only peak in the periodogram above 
the critical level of about 16 that corresponds to a false alarm 
probability of 0.01. Note that the false alarm probability does 
not take into account the flux density errors of the light curve. 

We also detect 2M0036H-18 with a flux density of 144 ± 
22 /iJy, similar to previous detections (©. No significant 
variability was detected in the light curve. 

3.2. VLBA observations 

TVLM513^6 is detected with the VLBA, the first such 
detection of any ultracool dwarf. To locate the source in 
a data set using the inner seven VLBA antennas we shifted 
the phase center of the uv data to the nominal VLA posi- 
tion and created a map covering 0.2" x 0.2", encompass- 
ing the VLA position uncertainty (Ict = 45 mas). We find 
a single robust source at the 8(t significance level, which 
persists with changes in the imaging parameters and source 
detection routines (AIPS SAD and JMFIT). The source is 
located 46 mas (i.e., \<j) from the nominal VLA posi- 
tion. We therefore identify this source as the VLBA coun- 
terpart of TVLM513^6 and shift the phase center to this 
position for subsequent analysis. The source is located 
at a position of RA=15'>01 "08. 1647696'^ (±0.0000077 s), 
Dec=+22°50'01. 63968" (±0.00011"), where the errors are 
statistical only (Figure[3]). 

With a synthesized beam size of 2.5 x 2.3 mas (PA = 26°) 
the source appears to be unresolved; it has a peak bright- 
ness of 420 ±48 /iJy beam"' and an integrated flux density of 
360 ± 75 ^mJy (the VLA and VLBA fluxes overlap at the 2a 
level). At the distance of TVLM5 13-46, the beam size corre- 
sponds to a spatial resolution of about 0.026 AU, or about 50 
times larger than the expected radius of TVLM5 1 3-46. 

To follow up on this detection with an increased angular 
resolution, we added the outer VLBA antennas located on 
Mauna Kea and St. Croix. The resulting synthesized beam 
size is 1.7 X 0.9 mas (PA = 6°). At this resolution, the S/N 
ratio decreases, and the peak brightness of TVLM513^6 
is 230 ±47 /iJy beam"', and its integrated flux density of 
540 ±150 /iJy has a larger error bar since the source no longer 
appears to be well-described by a single unresolved Gaus- 
sian component. Instead, it appears asymmetric and resolved 
(Figure |4]i. This indicates a marginally resolved source on 
the longest baselines, but the significance is too low to war- 
rant a firm conclusion. The synthesized beam corresponds to 
a source size of about 20 stellar radii, larger than expected 
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for the magnetic field scale (e.g., iLinskv & GarvllT983l) . We 
therefore conclude that if TVLM5 13-46 is indeed resolved, 
then this likely indicates a binary system with two radio- 
emitting sources and a projected separation of about 1 mas. 
In this case, since radio emission has only been dete cted to 
date f rom objects with spectral type earlier than L3.5 dBerged 
|2006|) . the binary would likely be close to equal mass. 

We note that 2M0036+18 was not detected in the VLBA 
data due to insufficient sensitivity. 

4. IMPLICATIONS OF THE VLBA DETECTION 

Since TVLM5 13-46 is unresolved on a scale of about 2.5 
mas, we can place a firm lower limit on the brightness tem- 
perature of its quiescent emission of > 4 x 10^ K. This is not 
sufficiently high to directly rule out thermal emission (e.g., 
IPulk 1985 ). To constrain the brightness temperature to values 
above 10^ K, indicative of non-thermal emission, the angular 
resolution will have to be < 0.4 mas. We note that other con- 
siderations (expected size of the magnetic field, radio spectra l 
index) point to non-thermal emission (e.g. JOsten et al.P 2Q06). 

The sub-milliarcsecond astrometry of TVLM5 13-46 that 
is available thanks to its detectability with the VLBA under- 
scores the potential for a VLBI astrometric companion search. 
Since TVLM5 13-46 is already at the stellar/sub-stellar tran- 
sition, such a compa nion woul d be a sub-stellar object (brown 
dwarf or planet). Bower et"aD ( 2009) have recently used such 
VLBA astrometry to search for companions to nearby early- 
M dwarfs. They detected five M dwarfs with spectral types 
M1-M5 in one to three epochs with the VLBA down to a 5(7 
sensitivity limit of about 0.5 mJy. Their astrometric residuals 
compared to optical measurements of the parallax and proper 
motion were about 0.2 mas, similar to the astrometric accu- 
racy we achieved here for TVLM513^6. 

The signal-to-noise ratio of our detection only allows us to 
search for an astrometric shift during the 7-hour observation 
when separately imaging the first and second halves of the 
data set. Using data from the inner 7 VLBA antennas we 
find a positional difference of 0.47 ±0.37 mas, consistent with 
no significant shift. We note that the expected maximum as- 
trometric shift (equal-mass companion) for an orbit of a few 
hours is about 0.2 mas, indicating that VLBA observations 
alone (with la ^ 0.2 mas) cannot detect the presence of such 
a companion at the brightness level of TVLM 513-64. How- 
ever, more sensitive combined VLBAh-GBT observations may 
provide the required signal-to-noise ratio to probe the pres- 
ence of an equal-mass companion in a few-hour orbit (corre- 
sponding to a separation of ~ few stellar radii). Moreover, if 
TVLM513^6 is indeed a roughly equal mass binary with a 
projected separation of about 1 mas (> 1 d orbit) we expect 
an easily detectable astrometric signature of ^ 1 mas on a 1 d 
timescale. 

Equally important, a planetary-mass companion to 
TVLM5 13-46 can be detected on longer timescales, since a 
wider orbit leads to a larger astrometric signal. In Figure|5]we 
plot the maximum reflex motion induced by a companion or- 
biting TVLM5 13-46 with a range of masses (1 - 40 My) and 
orbital periods (0.1 d to 1 yr). We find that at 3(t « 0.6 mas 
for VLBA observations similar to the one presented here, we 
can detect a 10 My companion with an orbital period of > 15 
d (> 0.06 AU), or a 2 My companion with an orbital period of 
> 0.5 yr (> 0.3 AU). Using the combined VLBAh-GBT, with 



an expected astrometric uncertainty of 3(7 « 0.3 mas (due to 
a higher signal-to-noise ratio and a higher fraction of longer 
baselines), we will be able to detect a 1 My planetary com- 
panion in a > 0.5 yr orbit. We note that these limits are ap- 
propriate for the case where TVLM513^6 is a single star If 
the object is an equal mass binary, as may be indicated by the 
results of direct imaging, the search for planetary-mass com- 
panions will be more challenging. 

5. CONCLUSIONS AND FUTURE PROSPECTS 

Using high angular resolution radio VLBI observations, we 
clearly detect the ultracool dwarf TVLM513^6 on an un- 
precedented small scale of about 2.5 mas or about 50 stel- 
lar radii. The source is unresolved on this scale, but it may 
be marginally resolved on a ^ 1 mas scale when using the 
longest VLBA baselines, although the S/N ratio is low. Given 
the corresponding physical scale of about 20 stellar radii, this 
may point to a binary system with two radio-emitting sources 
in a > 1 d orbit. 

Beyond the possibiUty of a spatially resolved binary, a 
search for an astrometric shift in the position of the source 
due to the gravitational influence of a close companion yields 
a 3a limit of < 0.6 mas, about 3 times larger than the ex- 
pected shift for an equal-mass companion in a few-hour or- 
bit. Since we lost the phased- VLA and GBT data, we expect 
that future VLBI observations will provide the requisite as- 
trometric accuracy to search for an equal-mass companion in 
such a tight orbit. These future observations will also allow us 
to determine whether the marginally-resolved emission from 
TVLM5 1 3^6 is indeed due to an equal mass companion with 
a projected separation of about 20 stellar radii. If this is in- 
deed the case, the expected maximum reflex motion is about 
1 mas, easily detectable with VLBAh-GBT observations that 
span several days. 

Equally important, the fact that TVLM513^6 is detected 
with the VLBA opens the possibility to astrometrically search 
for a planetary-mass companion with a > 10 d orbit. Mon- 
itoring of TVLM5 13-46 with the VLBA or VLBA+GBT on 
timescales of days to months to years will allow us to probe 
the existence of companions with masses of^ 1-lOMy and 
with orbits of ^ 15 d to ^ 1 yr. These observations may lead 
to the first detection of an extrasolar planet via radio astrome- 
try. 

Thus, the use of VLBI observations opens up a wide com- 
panion parameter space for TVLM513^6, with detection via 
direct imaging for a roughly equal-mass radio-emitting com- 
panion down to a scale of ^ 20 stellar radii, and detection via 
reflex motion on scales as small as a few stellar radii for an 
equal mass companion, and scales of > 0.1 AU for a plane- 
tary mass companion. This first VLBA detection also paves 
the way for similar observations of known and future radio- 
emitting ultracool dwarfs in search of companions. 
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Fig. 1 . — VLA 8.5 GHz light curve of TVLM513^6 with a time resolution of 10 s. The upper panel shows the total intensity (Stokes /), and the lower panel 
shows the circularly polarized flux (Stokes V); negative values indicate left-handed circular polarization. For better visibility, all points that have error bars more 
than three times larger than the average error were discarded. 
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P = 99% : 



Fig. 2. — Lomb-Scargle periodogram for the VLA 8.5 GHz light curve of TVLM513-46. The horizontal hne denotes the power level corresponding to a false 
alarm probabUity of 0.01. We find a significant peak at a period of 58.8 min. 



VLBA Observations of Ulti-acool Dwarfs 



7 



22 50 01.648 



o 
o 
o 

CM 



Z 

g 

z 
_i 
o 

a 




01.634 



01.632 



01.630 



15 01 08.1654 08.1652 08.1650 08.1648 08.1646 08.1644 08.1642 

RIGHT ASCENSION (J2000) 

Fig. 3. — VLBA detection of TVLM513^6 using the inner seven antennas. Contovir lines are -2,2,2 v'2,4,4v'2,8 a, where <t = 48 /nJy; the synthesized beam 
size is indicated in the lower left comer. 
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Fig. 4. — VLBA detection of TVLM513^6 using the entire VLBA (minus Hn). Contour lines are -2,2,2\/2,4,4v^ a, where a = 50 /iJy; the synthesized beam 
size is indicated in the lower left comer The source appears to be marginally resolved and asymmetric at this resolution, possibly indicating that TVLM513^6 
is a roughly equal-mass binary system. 
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Fig. 5. — Predicted maximum reflex motion for companions with masses of 80. 40, 20, iO, 5, 2, and 1 M; (left to right), as a function of the corresponding 
orbital period, assuming that TVLM5 1 3^6 is a single star. While a single observation with the VLB A+GBT has the potential to directly detect a radio-emitting 
(and hence roughly equal mass) companion down to a scale of ~ 1 mas (corresponding to a ~ 1 d orbit), observations spaced over a year can detect a sub-stellar 
companion down to a few Jupiter masses via reflex motion. If TVLM5 13-46 is indeed an equal-mass resolved binary on a milliarcsecond scale, then we should 
be able to further constrain the system with the detection of astrometric reflex motion, expected to be at a level of about 1 mas. 



